Coulomb scattering with remote continuum states in quantum dot devices 
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Electron capture and emission by Coulomb scattering in self- assembled quantum dot (QD) devices 
is studied theoretically. While the dependence of the Coulomb scattering (Auger) rates on the local 
wetting layer electron density has been a topic of intense research, we put special interest on the 
remote scattering between QD electrons and continuum electrons originating from a quantum well, 
doped bulk layers or metal contacts. Numerical effort is made to include all microscopic transitions 
between the Fermi distributed continuum states. The remote Coulomb scattering is investigated as 
a function of the electron density, the distance from the QDs and the temperature. Our results are 
compared with experimental observations, considering lifetime limitations in QD memory structures 
as well as the electron emission in pn-diodes. 

PACS numbers: 73.21.La,73.63.Kv,85.35.Be 
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I. INTRODUCTION 



The controversial discussions about the nature of the 
electron dynamics in quantum dots m (QDs) show that 
this topic is still not fully understood. Different electron- 
phonon coupling regimes leading to a multi-phonon pro- 
cess j^|3i,j^ or a polaron decay 5, 6] have been proposed. 
Especially for a multi-phonon process a so-called bottle- 
neck effect is expected due to the discrete QD density 
of states. 

On the other hand it has been shown theoretically 
IMI^^I as well as experimentally 11, 12] that electron 
capture and relaxation by Coulomb scattering (Auger ef- 
fect) with electrons in the wetting layer can be very effi- 
cient for large enough densities. Furthermore, an Auger- 
type process can overcome the bottleneck due to the con- 
tinuous energy dispersion in the wetting layer. 

The depletion of the wetting layer density by an electric 
field (1^ has been shown to reduce the relaxation rates. 
Similarly, the electron density in the vicinity of the QDs 
is reduced in a pn diode by the depletion layer within a 
DLTS experiment [Ij, [Ij, [Ig ■ Electrons which are can- 
didates for a Coulomb scattering process are not located 
at the QDs anymore, but at some distance. Therefore we 
are interested in Coulomb scattering with remote contin- 
uum electrons of a quantum well and bulk contact regions 
as a function of the distance from the QDs. 

We propose a theory to calculate the transition proba- 
bilities between two QD and two continuum states. De- 
bye screening effects are included according to the three- 
dimensional Poisson equation. From the two-particle 
scattering probability we calculate the total scattering 
probability of the electron relaxation within a QD (Auger 
process) and the excitation from a lower state into an ex- 
cited state (impact ionization or impact excitation) by a 
numerical integration over the Fermi distributed contin- 
uum states. 



II. THEORY 
A. General 

The band diagram of a QD and a remote quantum well 
and a bulk region defined by the doping profile is sketched 
in Fig. □ We assume that the 2DEG and the 3DEG 
are confined by infinitely high barriers. The transition 
probability per unit time between the QD states i and 
J, and the continuum states k and k' is calculated by 
Fermi's Golden Rule 
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where Ei^i and i^gn are the initial and final energies. The 
matrix element is given by 

(2) 

where tpijjl^, ipk and ipk' symbolize the respective wave 
functions ■ U represents the effective Coulomb poten- 
tial, which can be written in a two-dimensional Fourier 
representation as 



q' 



with 
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Here, and eg are the relative and absolute permit- 
tivities, e > is the elementary charge, A is the nor- 
malization area, z, z' are the coordinates in growth di- 
rection, r,!, r' are the in-plane coordinates, and q' are 
two-dimensional wave vectors with \q'\ = q' . The di- 
electric function accounts for static screening effects 
of free electron densities. The evaluation of the matrix 
elements (O depends on the QD wave functions, which 
themselves depend crucially on the QD shape. In order 
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FIG. 1: Sketched band diagram showing the remote Coulomb 
scattering with two- and three-dimensional continuum elec- 
trons 



to estimate the integral, we note that the wave functions 
tpi{r) and 4'jiT) are orthogonal and thus the lowest order 
contribution results from the dipole moment — ea of the 
effective charge density —etp*{r)'ipj{r). We approximate 
this dipole moment by two localized charges of opposite 
sign and thus obtain 

dVV;(r'){/(r, r')^.(r') « C/(r, |) - U{r, -|) , (5) 

where the center position of the QD is chosen as r' = 0. 
The total transition probability per unit time from a QD 
state i into a QD state j is then obtained from an inte- 
gration over the continuum states 
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where and f^i are the Fermi distribution functions 
of the initial and final continuum states. For a two- 
dimensional electron gas (2DEG) we substitute n = 2 
and n — A, and n — 3 and H. = AL for a three- 
dimensional electron gas (3DEG), where L is a (large) 
normalization length. 



B. Two-dimensional electron gas 

In order to find an expression for the dielectric function 
Ce for a two-dimensional electron gas positioned at 2 = d, 
we consider an elementary charge positioned at z = 0, 
corresponding to an external potential with the Fourier 
components 
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For the induced potential Ui^d originating from the in- 
duced electron density Sn2B (q) in the quantum well layer 



the Fourier components of the three-dimensional Poisson 
equation read as 



-g^C/ind(g, z) + d^Uindiq, z) = - 
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with the solution 
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Within the Debye screening approach, Sn2B{q) is given 
by 



dn 



Sn2Biq) = [Uindiq, d) + U,^t{q, d)] . (10) 

Assuming a single 2DEG quantization energy E2B, 
is given by 
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where Ep is the Fermi energy, m* is the effective mass, 
T is the temperature and ks is Boltzmann's constant, 
leading to the dielectric function for a 2DEG 
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with the inverse screening length 
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Inserting the two-dimensional Fourier expansion of the 
Coulomb potential into Eq. JSJ and using plane waves 
for the 2DEG states, we find 
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ttz and Oil are the z- and parallel components of a, k and 
k' are the two-dimensional wave vectors of the initial and 
final 2DEG states. Integration in real space gives 
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with the momentum transfer q = \q\ = |fe — 



C. Three-dimensional electron gas 

Coulomb scattering with semi-confined bulk electrons 
located at 2; > c? is considered next. The Fourier compo- 
nents of the three-dimensional Poisson equation for this 
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problem are 

-q^Umdiq, z) + 9f C/i„d(g, z) = - 



(16) 



where the induced bulk electron density is given by 
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With the definition of the Debye screening wave vector 
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and the expression for the bulk electron density 
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the Debye screening wave vector is given by 



( Ep — Ec \ 
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Here, denotes the Fermi integral of order s. Solving 
Eqn. lUni) leads to 
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Here, we have used that Uind is continuously differen- 
tiable at z = d. Analogously to Eqn. (|12|l we find the 
dielectric function for the inhomogeneous bulk electron 
distribution for z > d as 
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Using continuum wave functions, which are semi-confined 
in the z-direction and plane waves in the radial directions 
then leads to 
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where and fc^ denote the z-components of the wave 
vectors. Integration over real space results in 

fSDEG 1 2 _ 2e''e-29d [cosh((7az) - cofi{qa„)] 
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FIG. 2: Auger and impact excitation probabilities as a func- 
tion of the 2DEG electron density at d —0 nm (e.g., for a 
wetting layer surrounding the QDs) for different energy differ- 
ences AE and temperatures T. In addition the correspond- 
ing result without screening is shown for the Auger rate at 
AE =80 meV and 77 K (dash-double-dotted line). 



with qz = k'^ — kz and pz — k'^ + kz- Analogously to the 
2DEG, the transition probability is obtained by integra- 
tion over occupied continuum states, which is performed 
numerically. 



III. RESULTS 

A. Remote two-dimensional electron gas 

According to the geometry and size of InAs QDs which 
have a typical base length of about 10 nm and a strong 
quantization in z-direction, we use |a|||=5 nm and az=0. 
For small dipole moments a Taylor expansion in Eqn. 
(|15|l and (|24|l leads to a quadratic dependence on the 
transition dipole moments, and the transition probabili- 
ties presented here can be rescaled for specifically shaped 
QDs. 

First we performed calculations of the Auger and im- 
pact excitation probabilities per unit time as a function 
of the electron density per unit area in the wetting layer, 
using d = nm and m^^^^^ = 0.023 mg for different QD 
level spacings AE — \Ej — Ei\. Fig. [21 shows the result- 
ing probabilities for AE =30 meV and 80 meV for T=77 
K and 300 K. The Auger probabilities hardly depend on 
AE and the temperature. For densities less than 10^° 
cm~2 the Auger probability is a linear function of the 
density. From a fit Wi-,j = T'^gn2B for low densities 
we obtain an Auger relaxation coefficient of T^^^ =60 
cm2s^^. This value is about one order of magnitude 
larger than the corresponding value calculated in Ref. 
(1 cm2s~^). We explain this deviation by the different 
values for the transition dipole moment and the effective 
mass used in this work. Note that the Auger capture co- 
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FIG. 3: Remote Auger probabilities as a function of the dis- 
tance of the 2DEG from the QDs for n2D = lO^'^ cm~^, 
A£ =30 meV, AE =50 meV and AE =80 meV at 77 K 

efficients evaluated in Ref. (2x 10~^^ cm'^s""'^) are also 
much smaller than those observed experimentally in Ref. 
[Tlj l (10^^ cm'^s^^). For large densities, i.e., when the 
2DEG becomes degenerate, the density dependence be- 
comes nonlinear, and the probabilities are even reduced. 

Fig. El also shows that the impact excitation probabili- 
ties strongly depend on AE and T, and the ratio between 
the impact excitation and Auger probabilities is given 
by the Boltzmann factor exp{— AE / (kBT)) . The Auger 
probability for AE =80 meV and T=77 K but neglecting 
screening effects is also plotted in Fig. |21 demonstrating 
that screening effects play only a role for densities typi- 
cally larger than 10^° cm~^. 

Fig. O shows the remote Auger probabilities as a func- 
tion of the distance d between the QDs and the 2DEG, 
for = 10" cm-2, m^^^s = 0-06 toq, T = 77 K and 
AE =30 meV, AE =50 meV and AE =80 mcV. The 
Auger probabilities decrease exponentially with d, de- 
pending on the energy transfer. Therefore the Coulomb 
scattering of this type has a local character. As a conse- 
quence, a 2DEG at d =200 nm as in the memory struc- 
ture of Ref. 01 does not contribute to the lifetime limi- 
tation of charged QDs. 



B. Remote three-dimensional electron gas 

Next, we calculated the Coulomb scattering probabili- 
ties as a function of the bulk density in the vicinity of the 
QDs. The result is plotted for AE =30 meV (solid line), 
AE =50 meV (dashed line) and AE =80 meV (dotted 
line) in Fig. Hfor T=77 K, d =0 nm and TOg^As = 0.06 
mo- Again, for low densities the probabilities are linear 
in the density. But in comparison to the 2DEG, the fitted 
Auger coefficients = Wi-tj/n^Y) depend on the QD 
level spacings. Here we obtain T^^g=4.7x 10^^ cm'^s^^, 
=1.9x10-6 cm^s-i and T^^l.OxlQ-^ cm^s~^ for 



FIG. 4: Auger probabilities as a function of the three- 
dimensional electron density for different energy level sepa- 
rations AE at 77K and with d =0 nm. In addition the corre- 
sponding result without screening is shown for the Auger rate 
at AE =80 meV (dash-dotted line). 
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FIG. 5: The Auger probabilities for a three-dimensional elec- 
tron gas with nsD =10^'' cm"'^, m* =0.06 mo and the impact 
excitation probabilities for tisd =10^^ cm^'', m* = mo as a 
function of d for different level spacings AE at 77 K. Inset: 
dependence of the Auger probability for nsD =10^^ cm"^, 
m* =0.06 mo and AE = 80 meV on the transition dipole 
moment for d=100 nm: a^, 7^ 0, = (solid line); = 0, 
flz 7^ (dashed line) and ax=az 7^ (dash-dotted line) 



AE =30 meV, AE =50 meV and AE =80 meV, respec- 
tively. For high electron densities we observe a deviation 
from this linear behavior, since pairs of occupied and 
empty continuum states can only be found near the Fermi 
level when the electron gas becomes highly degenerate. 
The resulting Auger probability neglecting screening ef- 
fects for AE =80 meV is plotted as a dash-dotted line in 
Fig. ^ demonstrating that screening is only important 
for densities above 10^^ cm"^. 

The Auger probability as a function of the distance 
from the QDs is investigated next. The resulting curves 
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FIG. 6: Arrhenius plot of the inverse impact excitation prob- 
ability with d=90 nm, tisd ^SxlO^*^ cm"^ for AS =30 meV 
(solid line), AS =50 meV dashed line, A^: =80 meV (dotted 
line); dash-dotted lines: functions proportional to the Boltz- 
mann factor C exp(Aii'/(fcsT)) where C is chosen such that 
lines intersect at 1000/r=10 K"^ 

for nsD =10^^ cni-3, m^^^s =0-06 mo, AS =30 meV 
(solid line), A£; =50 meV dashed line, A£; =80 meV 
(dotted line) at 77 K are shown in Fig. [S] Compared 
to the remote 2DEG scattering, the dependence on d is 
weaker, since momentum in the 3DEG can additionally 
be changed in the third dimension. An interesting point 
is how a metal contact with its high electron density lim- 
its the lifetime of charged QDs. Therefore, the impact 
excitation probabilities are also displayed in Fig. [S] for 
n3D =10^'^ cm"'^, which is a typical electron density in 
metal contacts, and m* =mo. Depending on the energy 
spacing Aii^, we obtain transition times of the order of an 
hour for AE =80 meV and d =400 nm. This goes well 
with the observation of bistable behavior in QD struc- 
tures on this time scale [11 113. The inset of Fig. 
shows the dependence of the Auger probability on the 
transition dipole moment — ea, which can be well ap- 
proximated by a quadratic law resulting from the Taylor 
expansion of Eq. H23I) . The Taylor expansion also ex- 
plains that the transition probabilities are enlarged by 
a factor of 2 for vertical dipole moments (in growth di- 
rection) compared to the lateral dipole moments studied 
throughout this work. We conclude that the direction of 
the dipole moment does not affect the order of magnitude 
of the transition rates. 

Finally we consider QDs embedded in a pn diode, 
where the QDs are positioned within 3x10^^ n-doped 



GaAs layers. Such a device was investigated in DLTS 
experiments, and a level spacing of Ai? =82 meV with 
an emission time of 62 ms at T=40 K has been observed 
T^. The emission from the QD ground state level is 
dominated by an excitation into the excited state, from 
which the electron escapes into the continuum by a fast 
tunneling process. Taking into account the Boltzmann 
factor between excitation and relaxation processes we es- 
timate a timescale for the relaxation process from the 
excited state into the ground state of about 3 ps. From 
Fig. [21we conclude that for an Auger process on that fast 
time scale a wetting layer density of lO"'^" cm~^ would be 
necessary. However, the wetting layer is not occupied 
due to the electric field in the depletion layer. From the 
parabolic band bending and a QD ground state energy 
of 200 meV below the conduction band edge we estimate 
that the depletion layer extends to d =90 nm from the 
QDs. For these parameters, we display the temperature 
dependence as an Arrhenius plot as shown in Fig. I^lf^or 
/^E =30 meV (solid hne), /^E =50 meV dashed line, 
/S.E =80 meV (dotted line). The Botzmann factor (func- 
tions proportional to exp(Ai?/(fcBT)) chosen such that 
the lines intersect at 1000/r=10 K~^) is also displayed 
as dash-dotted line for each energy spacing, showing that 
the slope of the Arrhenius plot of the temperature depen- 
dence is mainly given by the energy spacing. 

Nevertheless, we obtain an impact excitation time of 
about 10"* s for T =40 K and AS =80 meV, which is 
orders of magnitudes slower than observed in the DLTS 
experiment. Only by reducing the level spacing below 
AE =50 meV impact excitation times in the range of 
tens of ms can be obtained at this temperature. This 
indicates that in addition phonon mechanisms have to 
be taken into account. 



IV. CONCLUSION 

In conclusion, we have evaluated the Auger effect and 
impact excitation (Coulomb scattering) in QD devices 
with remote continuum states. Three-dimensional De- 
bye screening has been included, which is important for 
high electron densities. Our results show that, depend- 
ing on the detailed device and QD geometries, scattering 
with remote continuum carriers can be of importance for 
the electron kinetics if the local wetting layer density is 
reduced. 
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